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ABSTRACT 

Peculiar velocities induce apparent line of sight displacements of galaxies in redshift 
space, distorting the pattern of clustering. On large scales, the amplitude of the distor¬ 
tion yields a measure of the dimensionless linear growth rate (3 « /b of fluctuations, 

where if the cosmological density and b is the bias factor. To make the maximum 
statistical use of the data in a wide-angle redshift survey, and for the greatest accu¬ 
racy, the spherical character of the distortion needs to be treated properly, rather than 
in the simpler plane-parallel approximation. In the linear regime, the redshift space 
correlation function is described by a spherical distortion operator acting on the true 
correlation function. It is pointed out here that there exists an operator, which is essen¬ 
tially the logarithmic derivative with respect to pair separation, which both commutes 
with the spherical distortion operator, and at the same time defines a characteristic 
scale of separation. The correlation function can be expanded in eigenfunctions of 
this operator, and these eigenfunctions are eigenfunctions of the distortion operator. 
Ratios of the observed amplitudes of the eigenfunctions yield measures of the linear 
growth rate /3 in a manner independent of the shape of the correlation function. More 
generally, the logarithmic derivative d/d\nr with respect to depth r, along with the 
square and component of the angular momentum operator, forms a complete 
set of commuting operators for the spherical distortion operator acting on the density. 
The eigenfunctions of this complete set of operators are spherical waves about the 
observer, with radial part lying in logarithmic real or Fourier space. 
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1 INTRODUCTION 

Peculiar velocities induce apparent line of sight displacements of galaxies in redshift space, distorting the pattern of clustering. 
In an influential paper, Kaiser (1987) pointed out that, in the large-scale linear regime, and in the plane-parallel approximation, 
the distortion takes a particularly simple form in Fourier space. He showed that a wave of amplitude 5{k) appears amplified 
in redshift space by a factor 1 -I- (in this paper a superscript s is used to distinguish the redshift space overdensity 6“ 
from the unredshifted overdensity S, and likewise the redshift space correlation function and power spectrum P'* from their 
unredshifted counterparts ^ and P): 

M(fe) = (l + d/r")5(fc) . (1) 

Here n = z.k is the cosine of the angle between the wavevector k and the line of sight z, and /3 is the growth rate of growing 
modes in linear theory, the dimensionless quantity which solves the linearized continuity equation V.n -I-/35 = 0 in units where 
the Hubble constant is one. In unbiased standard pressureless Friedmann cosmology, the linear growth rate /3 depends on the 
cosmological density parameter Q as (e.g. Peebles 1980, equation [14.8]) 

/3 « . (2) 


If the galaxy overdensity S is linearly biased by a factor b relative to the underlying matter density (5 m of the Universe, 
(5 = bSu, but velocities are unbiased, then the observed value of /? is modified to 




b 


(3) 


However, the assumption of linear bias lacks compelling justification, and equation (^) is really just an acknowledgment that 
the measurement of H may be biased if galaxies do not trace mass. In general, (3 depends on the adopted cosmological model, 
and is a function of H, the cosmological constant A, bias, and perhaps other quantities. 
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Kaiser concluded from equation (^) that the power spectrum of galaxies (the Fourier transform of the correlation function) 
appears in redshift space amplified by a factor (1 + 

PAk) = A + dnyPik). (4) 

Kaiser’s formulae and assume the plane-parallel approximation, where galaxies are taken to be sufficiently far 
away from the observer that the displacements induced by peculiar velocities are effectively parallel. To date, most studies of 
large-scale redshift space distortions have assumed the plane-parallel approximation (Lilje & Efstathiou 1989; McGill 1990; 
Loveday et al. 1992; Hamilton 1992, 1993a; Gramann, Cen & Bahcall 1993; Bromley 1994; Fry & Gaztafiaga 1994; Fisher et 
al. 1994a; Fisher 1995; Cole, Fisher & Weinberg 1994, 1995). In the plane-parallel approximation one is necessarily restricted 
to considering only pairs of galaxies separated by no more than some angle on the sky. For example, Hamilton (1993a) and 
Cole et al. (1995) restricted the opening angle to 50° as a reasonable compromise between statistical uncertainties and the 
error resulting from the plane-parallel approximation. Cole et al. (1994, fig. 8) show from simulations that at this opening 
angle the plane-parallel approximation causes /3 to be underestimated by about 5 per cent. 

Properly, however, the redshift displacements of galaxies are radial about the observer, not plane-parallel. A correct 
treatment of radial distortions, besides being more accurate for pairs of modest angular separation, would improve statistics 
by admitting pairs at large angular separations, which would be particularly helpful for measurements at the largest scales. 
In all sky surveys such as the IRAS redshift surveys, for example, one could effectively double (by going fore and aft instead 
of just fore) the effective scale over which redshift distortions can be measured. The need for greater accuracy will increase as 
redshift surveys grow larger, and statistical errors decrease. 

The first study of spherical redshift distortions was by Fisher, Scharf & Lahav (1994b) on the 1.2 Jy redshift survey. 
They expanded the density field in spherical harmonics, windowing the density in the radial direction with Gaussian windows 
at several depths. Heavens & Taylor (1995) improved on Fisher et al.’s procedure by expanding the radial density field in a 
complete set of spherical waves. Both sets of authors assumed a prior form of the (unredshifted) power spectrum in determining 
a maximum likelihood value of /3. Ballinger, Heavens & Taylor (1995) took Heavens & Taylor’s approach a step further by 
allowing the power spectrum to vary in six bins. One curious aspect of these studies is that they gave values consistently 
higher, P = 0.96jlQ ^g (Fisher et ah, 1.2 Jy), P = 1.1 ± 0.3 (Heavens & Taylor, 1.2 Jy), and P = 1.04 ± 0.3 (Ballinger et ah, 
1.2 Jy), than values obtained assuming the plane-parallel approximation, P = O.OOjljQ’j® (Hamilton 1993a, 2 Jy), P = 0.45j'lQ ^g 
(Fisher et al. 1994a, 1.2 Jy), and P = 0.52 ± 0.15 and P = 0.54 ± 0.3 (Cole et al. 1995, 1.2 Jy and QDOT). The origin of this 
discrepancy is not yet clear. 

The purpose of the present paper is to set forth a procedure for measuring P from spherical distortions in a manner 
independent of the form of the power spectrum. The basic idea is to represent the correlation function in eigenfunctions of 
the spherical distortion operator. 

It is instructive to see how this works in the plane-parallel case. In general, the linear redshift distortion equation (^) is 
an operator equation. In ordinary space, for example, it is 

= [1 J- /3(5/cIz)^V~^] 5{r) (5) 

where is the inverse Laplacian. A key advantage of Kaiser’s formulation is that Fourier modes are eigenfunctions of the 
(plane-parallel) distortion operator. This is intimately related to the circumstance that the iV operator, whose eigenfunctions 
define the Fourier modes, commutes with the distortion operator. 

The other key aspect of Kaiser’s formula (^ is that the unredshifted power spectrum P{k) is a function only of the 
absolute value k of the wavevector. This implies that ratios of the redshift power spectrum for waves with the same k but 
different angles to the line of sight (different /r) yield measures of P independent of the amplitude P{k) of the power spectrum. 
This argument identifies k^, or equivalently the operator — V^, as playing a special role, which is, roughly stated, that k^ 
‘defines a scale’ of separation. 

The first aim, then, of this paper is to identify an operator, analogous to k^ in the plane-parallel case, which both 
commutes with the spherical distortion operator, and at the same time ‘defines a scale’ of separations. This is done in Section 

3, the spherical distortion operator having been derived in Section 2. The procedure is carried a step further in Section 

4, which presents a complete set of commuting operators for the spherical distortion operator. The eigenfunctions of this 
complete set of commuting operators are logarithmic spherical waves. Section 5 discusses how to measure p. The conclusions 
are summarized in Section 6. 


2 THE SPHERICAL DISTORTION OPERATOR 

This section derives the operator equations © and ( p^ which relate the redshift density, hence correlation function, to the 
real density and correlation function in the linear regime. We assume the standard gravitational instability picture in the 
standard pressureless Friedmann cosmology (e.g. Peebles 1980). The derivation treats correctly the spherical character of the 
redshift space distortions about the observer, so structures may be nearby and may subtend a large angl e. It is nevertheless 
necessary to exclude a local region of the Universe about the Milky Way, both to avoid local bias (Section 2.1), and to ensure 
the linear requirement that peculiar velocities be small compared to the distance to the observed structures, equation (|li]). 
The derivation here is similar to that of Kaiser (1987, section 2). 
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In this subsection, the frame of reference is taken to be stationary, that is, the frame of reference of the Cosmic Microwave 
Background (CMB). Subsection 2.2 below discusses the transformation to the Local Group frame. 

Let s denote the observed redshift distance (in the CMB frame) to a galaxy, and r its true distance from the observer. 
The two distances differ by the line of sight peculiar velocity v of the observed galaxy: 


s = r + n (6) 

in units where the Hubble constant is unity. The number N'‘{s)A^s of galaxies observed in an interval d®s of redshift space in 
a redshift survey is related to the real space number density N{r) by number conservation: 

N‘’{s)Ah = N{r)A\ . (7) 

An unbiased estimate of the true galaxy overdensity 5(r) at position r is given by 1 + S{r) = N{r)/^{r), where $(r) is the 
selection function, the expected number density of galaxies at depth r given the selection criteria of the survey. As emphasized 
by Fisher et al. (1994b), the selection function <E>(r) in a flux-limited survey is a function of the true distance r, not of the 
redshift distance s. Thus one would be inclined to estimate the overdensity (5'*(s) in redshift space by l-|-5'’(s) = A'^'*(s)/>l>(r). 
However, this estimate of the overdensity 5'’(s) requires knowing not only the true selection function ^{r), but also the true 
distance r to the galaxy at redshift distance s. This can be done, but it requires a full reconstruction of the deredshifted 
density held, which is not the philosophy of the present paper. We assume here instead that one measures, by some standard 
technique (Binggeli, Sandage & Tammann 1988), a selection function 4>'’(s) in redshift space, and that the overdensity S‘‘{s) in 
redshift space is then dehned by 1 -|- 5'’(s) = N‘‘{s)/$‘‘{s). We distinguish here between the measured redshift space selection 
function 'l?'’(s) and the true selection function ^{r), though, as argued below, the two in fact agree to linear order. Thus the 
relation between the observed redshift space overdensity 5'’(s) and the true overdensity S{r) is 


'l>'’(s)[l -I- (5'’(s)]s^ds = ‘I>(r)[l -I- S{r)]r^dr . 

(8) 

With equation (^, equation (^ rearranges to 



( 9 ) 


which, with some small differences, is Kaiser’s (1987) equation (3.2). In linear theory, the peculiar velocity of the galaxy along 
the line of sight is 

v{r) = -(3^y~'^S{r) (10) 


where V ^ denotes the inverse Laplacian. Thus to linear order in the overdensity S{r), and with the additional condition that 
peculiar velocities v be much smaller than the distance r to the structures being observed. 


u < r , 
equation 

SSr) = 


reduces to [note in particular that ^“(s) = 5‘’{r) to linear order] 


^+d{^ + 


i{r)d 

rdr 


S{r) 


( 11 ) 


( 12 ) 


where a(r) is the logarithmic slope of times the redshift space selection function 'l‘'’(r) at depth r: 


a{r) 


d In (r) 
91nr 


(13) 


Equation (^) is an operator equation relating the redshift space overdensity 5“ to the real space overdensity 5. The quantity 
in square brackets in equation (|^) is the spherical distortion operator. It is straightforward to transform the equation 
into Fourier space (Zaroubi & Hoffman 1995), but, unlike the plane-parallel case considered by Kaiser (1987), the spherical 
distortion operator does not simplify to an eigenvalue in Fourier space, but remains an operator. Tegmark & Bromley (1995) 
obtain an expression for the Green’s function (i.e. the inverse) of the spherical distortion operator in equation (|l^), for the 
particular case a = 2. 

Deriving the linearized equation ( p^ ) from equation (^ involves the approximation <E>'’(r) = 4?(r), which is valid to linear 
order, as will now be argued. In a flux-limited survey, the selection function <l>(r) at depth r is the number density of galaxies 
luminous enough to be seen to depth r, which is the integrated luminosity function above a threshold. All modern methods of 
measuring the selection function (Binggeli et al. 1988) seek to eliminate dependence on density inhomogeneity, by assuming 
that the luminosity function is independent of density, and constructing the luminosity function by counting relative numbers 
of bright and faint galaxies in identical volumes. If bright and faint galaxies at any point share the same peculiar velocity, as 
seems probable in the linear regime, then the measurement of the selection function in redshift space should remain unbiased 
by inhomogeneity. A systematic difference between the real and redshift space selection functions will however result from the 
fact that on average peculiar velocities will cause galaxies to tend to ‘diffuse’ away from depths where the density per unit 
velocity r^<I>(r) is larger, towards regions, both shallower and deeper, where the density is smaller. This diffusion will perturb 
the shape of the luminosity function observed in redshift space from its true shape, biasing the measurement of 4?“(r). The 
diffusive flux of galaxies is proportional to the variance (Au^) of radial peculiar velocities at depth r, which is of second order 
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in perturbations, and hence vanishes to linear order. It follows that the redshift and real selection functions agree to linear 
order, = ‘I>(r), as claimed. It is to be noted that in general the difference between il>“(r) and <I>(r) is confined only to the 

monopole mode about the observer, and only to the longest wavelength modes, since the difference of the functions is slowly 
varying with depth r. Thus in any case there is little harm in approximating $'’(r) = <E>(r) for all but the fundamental modes. 
Fisher et al. (1994b) make the clever point that the observed galaxy density weighted by some window W(s) is 

W{s)N‘’{s)dds = W{s)N{r)<fr « IT(r) (l + N{r)d\ (14) 

\ r alnr / 

for sufficiently smooth windows W{s). Equation (^^, which is basically an integration by parts, allows the effect of redshift 
distortions to be cast in terms of the window IF(s), rather than the selection function 4>“(s) and the overdensity S‘’{s). 
However, our aim here is to decouple the measurement of the linear growth rate /3 from the shape of the correlation function. 
For this purpose it appears essential to work with the overdensity il®(^ deconvolved from the shape of the selection function, 
which compels us to work with equation ( [l2| ) rather than equation (Q). 

The redshift space correlation function follows from equation (h^by taking an ensemble average of products of densities 
at positions r\ and r 2 relative to the observer: 


5''(r-i2,ri,r2) = (5'*(ri)5''(r2)) 


1 + /3 


a(ri)d\ 
drf ndn J 



1 + P 


a(r2)9\ 
dr2 r2dr2 J 



i{ri2) 


(15) 


where ri 2 = \r\—r 2 \ is the pair separation. Equation (^^ is the basic equation which describes the redshift space correlation 
function ^'*(ri 2 ,ri,r 2 ) in the linear regime as the result of the spherical distortion operator acting on the true correlation 
function ^(ri 2 ). It is a consequence of isotropy about the observer that the redshift correlation function ^'’(ri 2 ,ri,r 2 ) is a 
function only of the lengths of the sides, not of the orientation, of the triangle dehned by the observer and the pair of galaxies 
observed. 


2.1 Local bias 

Equation (^|) is the ensemble average redshift correlation function observed by stationary observers at random positions in 
the Universe. But we are not at a random position: we reside on a galaxy, the Milky Way. 

Consider an ensemble of observers who sit on galaxies, measuring the correlation function. The probability that an 
observer on a galaxy at position 0 sees a pair of galaxies at points 1 and 2 is the probability of finding a triple of galaxies at 0, 
1, and 2, divided by the probability of finding the galaxy at 0. It follows that the apparent correlation function 5i2(apparent) 
observed from galaxies at 0 is 

Ci 2 (apparent) = ((1 + (5o)5i(l2) = C 12 + C 012 (16) 

which differs from the true correlation function ^12 by the three-point correlation function ^ 012 . The three-point function in 
the real Universe appears well approximated by the hierarchical model (012 = Q(€oi€i 2 -I- 502 C 12 +C 01 C 02 ) with Q « 1 (e.g. Fry 
& Gaztanaga 1994, table 8). To ensure that the three-point term in equation ( [l6| ) is small compared with the desired two-point 
term, <^012 ^ ^ 12 , requires at least that ^ 01 , ^02 ^ 1- Thus, on average, galaxy-bound observers who want an unbiased measure 
of the correlation function would be advised to exclude a local region around themselves which is at least a correlation length 
in radius. 

The importance of deleting the local region around the Milky Way needs emphasizing. If the local region is not deleted, 
then the spherical distortion equation (^^ will attempt to interpret the local overdensity as caused by peculiar redshifts, 
which is clearly wrong. The local region must be deleted in any case to ensure that peculiar velocities are small compared 
with the depth of the structures being observed, u <C r, equation (pd]). 


2.2 The peculiar velocity of the Milky Way 

Formula ( p^ ) for the redshifted correlation function is valid for randomly located observers in stationary frames of reference. 
The Milky Way however is not stationary, but moving. Happily, this motion is rather accurately known from the dipole 
anisotropy of the CMB (Kogut et al. 1993). The linear part of the Milky Way’s peculiar velocity is the streaming motion V 
of the Local Group (LG) (Yahil, Tammann & Sandage 1977) with respect to the CMB frame. 

In the LG frame, the peculiar redshift of a galaxy at true position r relative to the observer is 

s'"® = r + v — r.V . (17) 

Recapitulating the derivation of equation but starting from equation (|^) instead of (^, one concludes that, to linear 
order, the redshift space overdensity 5®'"®(r) observed in the LG frame differs from the redshift space overdensity ^“(r) in 
the GMB frame by a dipole term directed along the LG motion V: 

^sLG(^) _ . (Ig) 

r 

The fact that the dipole term changes sign when a{r) goes negative, that is, where the selection function 'l'(r) is steeper 
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than r is the origin of the ‘rocket effect’ noted by Kaiser (1987, section 2). It follows from equation ( |l^ that the redshift 
correlation function ri,r 2 ) in the LG frame is related to that 5 '’(ri 2 ,ri,r 2 ) in the CMB frame by 


r^°(n2,ri,r2) = r(ri2,ri,r2) + 


a{ri)V5 


iLG 


3ri 


(riV) ^2)V5^ 


LG 


3r2 


{r2V) a{ri)a{r2)V‘^ 


3rir2 


ri.r2 


(19) 


where 5i'^®(rV) is the dipole component of the redshift space overdensity in the LG frame at depth r in the direction of the LG 
motion V. That is, if the overdensity is expanded in spherical harmonics about the observer, = y~),^ (r)Yim(r), 

then (5f^''(r) = i® ^^e dipole component. As in the case of the density, the difference between the correlation 

functions in the two frames is pure dipole. 

There are basically three ways to deal with the finite peculiar velocity V of the Local Group with respect to the CMB 
frame, when measuring redshift distortions. The first is to work in the CMB frame using the known value of V. The second 
is to treat V as unknown, and to average the right hand side of equation (^^ over an ensemble of observers who measure 
structure in frames at rest with respect to the local streaming frame. The third way is to avoid the problem altogether by 
ignoring the dipole component of the redshift correlation function. 

Of these options the first, to work in the CMB frame using the known value of V, makes the most use of available 
information, and is also the simplest, which would seem to make it a clear winner. 

To linear order in the supposedly small quantity v/r, measuring in the CMB frame is the same as measuring in 
the LG frame and then using equation ( p^ to transform to the CMB frame. The two measures do however differ in higher 
order. Since for nearby structure the approximation v — r.V <C r is better than the approximation u <C r, it should be 
somewhat more accurate to measure in the LG frame and then to transform to the CMB frame using equation (p^). 


3 THE a/a In ri 2 1A OPERATOR 

The symbol A is shorthand for any two variables, such as (ri/ri 2 , r 2 /ri 2 ), which define the shape (angles) of the triangle 
formed by the observer (= the Milky Way) and the observed pair of galaxies at positions ri and r 2 relative to the observer. 


3.1 The constancy of a in the distortion operator 

The spherical distortion operator, the quantity in square brackets in equation (p^), contains a term which depends on the 
function a{r), equation (y_^, which is the logarithmic slope of times the selection function >l>(r) at depth r. It is a basic 
assumption of this paper that a(r) can be approximated by a constant (or more generally by a function only of triangle 
shape A), to ensure that the operator <?/91nri2|A considered in Section below, and also the logarithmic radial derivative 
operators considered in Section 4, do in fact commute with the spherical distortion operator as claimed. 

While a(r) is generally not constant in real redshift surveys, it is typically a slowly varying function of depth r. According 
to the uncertainty principle, modes can be located in space within no better than a wavelength, so the approximation of 
constant a might be expected to be valid for modes whose wavelengths are short compared with the scale over which a{r) 
varies. 

In a flux-limited redshift survey, a is typically approximately constant at moderate depths, a consequence of the power- 
law character of the luminosity function at faint fluxes. Since a decreases in importance for pairs at depths greater than their 
separation, the approximation of constant a should be valid at least at separations no larger than the depth at which a starts 
to deviate significantly from a constant. 


3.2 The d/d\nr\ 2 \t\ operator 

Let 9/91nri2|A denote the logarithmic derivative with respect to pair separation ri 2 , the shape (angles) of the triangle formed 
by the observer and the observed pair of galaxies being held fixed. 

The operator 9/91nri2jA possesses two key properties, discussed in the Introduction. The first is that it commutes 
with the spherical distortion operator (to the extent that a is constant or can be approximated as a function of A, Section 
3.1), which implies that eigenfunctions of d/d\nri 2 \A are also eigenfunctions of the spherical distortion operator. The second 
property is that the eigenfunctions of 9/91nri2|A form a complete set for the (unredshifted) correlation function ^(ri 2 ), which 
is a function of separation ri 2 alone. 

A complete set of orthogonal eigenfunctions of d/d\nri 2 \A is 

( 20 ) 

where aJi 2 takes all real values from negative to positive infinity, and 7 is some fixed real numb^ chosen in practice to secure 
convergence at ri 2 —> 0 and ri 2 ^ 00 . The choice of the units of the separation ri 2 in equation (M) is a matter of convenience. 
Let 7 i 2 denote the representation of the real space correlation function ^(^ 12 ) in the eigenfunctions 


i{ri2) = Co. 


McJi2 


7i 2 = (27r)' 


C(’'i2)r7'“^^dri2/ri2 


( 21 ) 
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In effect, is the Fourier transform of ^{ri2)rj2 with respect to logarithmic separation lnri2. If the index 7 + iuJi2 were 
taken to be real rather than complex, then would be the Mellin transform of ^(ri2). The reality of 5(?'i2) implies the 
condition 




( 22 ) 

Then the correlation 


Consider for example the case where the correlation function is an exact power law ^(ri2) = rH- unc i^uiiciauiun 

function in aii2-space is a Dirac delta-function, = ( 5 d(<ui 2 ), provided that 7 in the expansion (E^ is chosen equal to the 

actual power-law index 7 of the correlation function. 

More realistically, if the correlation function is flatter than rH a-t small scales and steeper than at larger scales, then 
the correlation function in a;i2-space should look like a ‘window’ about tJi2 = 0 of width ~ 1 /Alnr, where Alnr is roughly 
the logarithmic range of ri2 over which ^ approximates a power law r^- 

The representation ^'*(a;i2. A) in a;i2-space of the redshift correlation function is defined similarly to equation 


/ oo noo 

^'’(0112, A)r■'^■“‘' 1^0112 , ,t“(a)i 2 . A) = (27r)“M ^(^12, A)r 72 ^“‘'iMri 2 /ri 2 |A • 

- 00 0 

In aii2-space, the spherical distortion equation becomes 
C*('^i 25 A) = [1 -I- / 3 Ai(lji 2 , A) -|- / 3 ^A 2 (aii 2 , A)] 
where, with 


712 = 7 -I- ia;i2 , 
the coefficient of /3 is 

^^3 3(3-712) 

while the coefficient of 0^ is 


B2(A)-h 


( 3 - 712 ) 


Bi(A) 


( 23 ) 

( 24 ) 

( 25 ) 

( 26 ) 


A2(a72,A) = i-^^i^i32(A) + 


1-b' 


5 - 712 \ 2 - 712 

The shape-functions Si (A) in equations (| 
Bo (A) = 1 


8712 (2 -|- 712) 
35(3 -7i2)(5-712) 

Bi(A) _ / 


B4(A) - 


2(1-712) 


(3-7i2)(5-7i2) 


B3(A) 


B3(A) 


3 - 712 ^ 

and (M) are (with a constant function Bo to complete the set) 




ri 


r2 


rW2 ^ ^ ^ _r?2/^ ^ 2(1 - 7f2) 


rir2 


rir2 


B2(A) = -[P2(/il)+B2(/r2)] 

B3(A) = 1 — 102 

B4,{A) = i(35/r?/ri — 15/1? - 15/ri -I- 3) 


( 27 ) 

( 28 ) 

( 29 ) 

( 30 ) 

(31) 

( 32 ) 


where P2{0) = |/i^ — | is the quadrupole Legendre polynomial, and 712 = cos0i2, 71 = cos61 and 72 = cos02 are the cosines 
of the interior angles of the triangle formed by the observer and the observed pair of galaxies, as illustrated in Fig. 


712 


2,2 2 

ri -I- r2 - ri2 


Ml 


2.2 2 
ri2 +ri-r2 


M2 = 


2.2 2 
ri2 +r2-ri 


2rir2 ’ 2ri2ri ’ 2ri2r2 

The plane-parallel limit is attained when the distance to each galaxy of a pair is large compared with their separation, 
ri, r2 S> ri2. Here 71 —» —72 and 712 —> 1 , the functions Bi and B3 vanish, and the functions B2 and B4 go over to Legendre 
polynomials B2(A) —> s-nd B 0 A) P^in) in the cosine I71I « \fi2\ of the angle to the line of sight. Thus the 

spherical distortion coefficients Ai and A2 in equations (^) and (g^ go over to their plane-parallel limits (compare Hamilton 
1992 , equations [ 10 ]-[ 12 ]) 

2 4712 p 

3 3(3-712)^"^^^ 

1 4712 p ( \ : 8712(2-1-712) 

5 7(3-712) 35 ( 3 -7i2)(5-712) 

In the plane-parallel limit, the distortion decomposes naturally into a sum of mutually orthogonal parts, the harmonics Pi{n), 
but in the general case the shape-functions Bi{A), equations (| 2 ^)-(| 3 ^), are not orthogonal. 

In the opposite limit where the distance to one of the galaxies, say 1 , of a pair is much smaller than their separation. 


Ai 


A 2 


P4(m) • 


( 34 ) 


( 35 ) 
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Figure 1. Contour diagrams showing, on the left, the ‘quadrupole’ amplitude Ai —2/3, and, on the right, the ‘hexadecapole’ amplitude 
^2 — 3/7Ai+3/35, equations (26), (27) and (38), of the spherical distortion of the redshift correlation function for the case of a power-law 
correlation function ^ oc r -^2 with index 7 = 1.5. To demonstrate the influence of the terms depending on the shape of the selection 
function, i.e. on a, the top two contour diagrams are for a = 2, appropriate for a volume-limited sample, while the bottom two are for 
q; = 0, appropriate for a selection function $(r) oc r~^. The units are such that the unredshifted correlation function has unit amplitude. 
The geometry of the observer and the observed pair is illustrated by the stick diagram at the centre. In the diagrams, the observed pair, 
represented by two large dots joined by a thick vertical line, is flxed, while the position of the observer varies over the diagram. Thus 
observers at 6 or 12 o’clock observe pairs which are aligned along the line of sight, while observers at 3 or 9 o’clock observe pairs which 
are transverse to the line of sight. The kink in the contours at the outermost radius in each diagram is not an error; the final ring shows 
the plane-parallel limit, corresponding to infinite radius. In the ‘quadrupole’ distortion, on the left, the correlation function appears 
reduced for pairs along the line of sight, and enhanced for pairs transverse to the line of sight, which is the expected squashing effect. 
If one imagines the figures to be rotated about their vertical axes, then volume elements in the resulting three-dimensional figure are in 
correct proportion to the volume element dA = ■ridrir 2 dr 2 /r ^2 geometrical configuration. Greyscales are graduated at intervals 

of 0.1 in the interval [—2,2], with contours at intervals of 0.5. Outside [—2,2], contours are at intervals of 1, up to ±10. 


© 1995 RAS, MNRAS 000, 000-000 



8 A. J. S. Hamilton and M. Culhane 


ri <C ri 2 , then r 2 ^ ri 2 , ^12 
and © reduce to 


. a ri2Mi 

^ (3 - 7712 ) ri 


—/ii and 7^2 ^ 1, the Bi terms become large compared with the others, and equations (|26| 


(36) 


A2 


a 

5 - 7?12 



2 - rii2 


ri2Mi 

ri 


(37) 


Here the distortion is dominated by a dipole about the near galaxy, reflecting the streaming motion of galaxies past the 
stationary observer. The divergence of the distortion as ri —> 0 is an artefact which results from a violation of the presumption 
( p^ that V n. The divergence is removed in practice by deleting the local region. Note that the validity of the linear 
approximation requires not that the dipole term be small, but only the weaker constraint that ^(ri 2 )ri 2 /ri be small. 


3.3 Picture 


Fig. 1^ illustrates the spherical distortion of the redshift correlation function for the case of a pure power law ^ oc r ^2 with 
index 7 = 1.5. As mentioned above, ^ui -12 is ^ Dirac delta-function at 0112 = 0 for a pure power law, so 7712 = 7 = 1.5 in 
equations (^^ and ( [ 2 ^ for this case. Fig. ^ shows not the spherical distortion coefficients Ai and A 2 directly, but rather 
the combinations Ai — 2/3 and A 2 — SjTAi -|- 3/35 which go over respectively to pure quadrupole and pure hexadecapole 
distortions in the plane-parallel limit (cf. equations and j^). In terms of these coefficients, which we refer to loosely as 
‘quadrupole’ and ‘hexadecapole’ coefficients, the distortion equation (^^ is (the normalization of the coefficients here differs 
by factors of 3/4 and 35/8 from the normalization of the ‘quadrupole’ and ‘hexadecapole’ correlation functions H 2 and of 
equations Q & in Section |5.2|) 


r 





3, 3 

7 ^' + ^ 



(38) 


To demonstrate the influence of the terms depending on the shape of the selection function, i.e. on a, equation (p^. Fig. 
1 ^ shows two cases, the first with a = 2 , appropriate for a volume-limited sample, the second with a = 0 , appropriate for a 
selection function 4>(r) oc r~^. 

One of the striking aspects of Fig. is how large the distortions become for ‘nearby’ pairs, those whose distances are less 
than or comparable to their separation, r\, r 2 r \2 (note that such pairs need not be physically nearby if the separation 

ri 2 is large). The large nearby distortion is produced by the ‘dipole’ term Bi(A), equation (j2^. This dipole should not be 
confused with the dipole in the correlation function, whose general form is /ii 2 times some function of ri and r 2 , independent 
of separation ri 2 . 

The amplitude of the nearby dipole distortion depends on a. This raises the concern that approximating a by a constant 
[or more generally by a function Q!(A)] may be unreasonable, if in fact the amplitude of the distortion depends sensitively 
on a. However, the problem is not as bad as it seems, because the functions Bi for i ^ 1 are all orthogonal to B\ in the 
limit where one galaxy of a pair is much closer than their separation. That is, Bi for i 7 ^ 1 all go over to sums of a constant 
and a quadrupole in the limit ri/ri 2 or r 2 /ri 2 —> 0, so are orthogonal to the diverging dipole. Thus uncertainty caused by 
approximating a translates mainly into uncertainty about the amplitude of the B\ part of the distortion, with lesser impact 
on the other parts. 

It is to be noted that the spherical distortion coefficients Ai and A 2 , equations (^^ and (© , depend on a only through 
the shape-functions Bi and B 3 . As will be seen in Section this has the consequence that the spherical distortion separates 
naturally into a part (Hf and Hf, equations [Q and |^^) that depends on a, and a part (HJ, and H|, equations [f77|-[[79t) 
that has no explicit dependence on a. 


4 LOGARITHMIC SPHERICAL WAVES 

The trick in Section 3 of expanding the correlation function in eigenfunctions of a single operator fails in Fourier space. For 
example, one would naturally think to expand the (unredshifted) power spectrum P{k), which is the Fourier transform of 
the correlation function ^(ri 2 ), in eigenfunctions of an operator d/dlnk. The procedure fails because the power spectrum 
is correctly a function 5u{ki + k 2 )P{k\) of two arguments, not one. The Dirac delta-function Suiki -I- fe 2 ) arises from the 
statistical homogeneity, or translational invariance, of clustering in space, the Fourier modes 5{k) being by definition the 
eigenfunctions of the translation operator iV. In redshift space radial distortions destroy homogeneity, although statistical 
isotropy about the observer is preserved. Thus the delta-function ^(fci -I- fe 2 ) cannot be an eigenfunction of the spherical 
distortion operator, and it cannot be factored out. The situation here contrasts with the plane-parallel limit, where translation 
symmetry is preserved in redshift space, and the delta-function factors out of the distortion equation. 

So is there a basis of eigenfunctions for spherical distortions analogous to the Fourier modes for plane-parallel distortions? 
As just mentioned, Fourier modes themselves fail, precisely because spherical distortions destroy homogeneity. 

Now the great advantage of the power spectrum is that, for Gaussian fluctuations, the (unredshifted) modes 5{k) are 
independent. Their independence is a consequence of homogeneity coupled with the Gaussian assumption. Homogeneity 
implies that modes with different wavevectors k are uncorrelated, {S{ki)5*{k2)) — 0 for fci 7 ^ k 2 , while Gaussianity means 
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that all higher order correlations vanish. Linear fluctuations may well be Gaussian, a consequence of an alliance between the 
Central Limit Theorem and physical smoothing processes in the early Universe. 

Spherical distortions destroy homogeneity, but they preserve isotropy about the observer. Thus the advantage of mode 
independence is preserved to the greatest extent if the density is expanded in eigenfunctions of the generator of rotations, 
which is the angular momentum operator L (e.g. Landau & Lifshitz 1958). These eigenfunctions are spherical harmonics, 
and symmetry about the observer implies that spherical harmonic modes with different Im will be independent for 

Gaussian fluctuations. This basic advantage of the spherical harmonic modes was exploited by Fisher et al. (1994b), and by 
Heavens & Taylor (1995). 

This leaves the problem of finding radial eigenmodes for the spherical disto rtion operator. Here one notices that, at least 
to the extent that the slowly varying function a{r) is constant (see Section 3.1), the spherical distortion operator in equation 
(@) is scale-free. A consequence of this is that the logarithmic derivative d/dlnr with respect to depth r (or alternatively 
the logarithmic derivative d/d\nk with respect to radial wavevector k) commutes with the spherical distortion operator, and 
its eigenfunctions therefore provide a basis of rad ial e igenmodes. 


We proceed with this idea further in Section 4.2, but first it is necessary to write down some standard definitions. 


4.1 Fourier space definitions 

Let 5“(fc) denote the Fourier transform of the overdensity: 


5Sk) = (Stt) 


-3/2 


5“(r)e''“’-dV 


&Sr) = (27r) 


-3/2 




(39) 


The (unredshifted) power spectrum P(fc) is by definition the Fourier transform of the correlation function ^(ri 2 ), with 
the conventional normalization 


P{k) = / ^(ri2)e‘'' ’’^^d®ri2 


^(^■ 12 ) = ((5(ri)(5(r2)) = (27r)' 


P(fc)e“‘'' ’'i=d®fe . 


(40) 


In redshift space, it is necessary to retain the dependence of the correlation function ^“(ri, r 2 ) and its Fourier transform 
the power spectrum ^‘‘{ki,k 2 ) on both their arguments: 


f (fci,fc2) = {5%ki)S%k2)) = (2^)“" j 


(41) 

^'’(r’i,r2) = (5'’(ri)5'’(r’2)) = {2-k)~^ J 

f (fci,fc2)e"‘('“i "i+''" ""^d%d^fc2 . 

(42) 

The unredshifted power spectrum ^(fci 

,fc 2 ) is related to the ‘reduced’ power spectrum P{k), equation (^^, by 


|(fel,fe2) =<5D(fel-hfe2)P(fel) . 


(43) 


4.2 A complete set of commuting operators 

In Section ^ it was shown that, for linear fluctuations, the overdensity 3‘‘{r) in redshift space is described by a spherical 
distortion operator acting on the real overdensity 5{r), equation (^^. A complete set of commuting operators (to the extent 
that the quantity a is constant. Section 3.1) for the spherical distortion operator is 

d 


dlnr 


d 

dink 


-3 


Lz 


(44) 


Here d/dlnr is the logarithmic derivative with respect to depth r, which is the same, up to a change of sign and a constant, 
as the logarithmic derivative d/dlnk with respect to radial wavevector k in Fourier space. The operators and Lz are the 
square and 2 -component (along some arbitrary axis) of the angular momentum operator L — ir x d/dr = ifc x d/dk, which 
is the same operator in real and Fourier space. The eigenfunctions of d/dlnr or d/dlnk are radial waves in logarithmic 
depth r or wavevector k, while the eigenfunctions of and Lz are the usual orthonormal spherical harmonics Yim- Thus the 
eigenfunctions of the commuting set (^^ are spherical waves with radial parts in logarithmic real or Fourier space. 

Equations (^)-(^) below are valid for both redshifted and unredshifted correlation functions. Let denote the 
representation of the overdensity S‘‘{r) as spherical waves in logarithmic real space: 


SSr) = (27r) 


- 1/2 




-7/2-i“y 

ri' In 


,(r)da 


= (2^) 


- 1/2 


S‘{r)r 


— 3 + 7/2 + iu; 


yi*m{r)d\ 


(45) 


and let 5^;^ denote the alternative representation of the overdensity as spherical waves in logarithmic Fourier space 

SSk) = i2n)-y^ f Y,6‘^,znk-^+y^+'^Y^{k)du; , = f 6‘{k)k-y^-'-YUk)d^k . 

—00 J 


(46) 


© 1995 RAS, MNRAS 000, 000-000 







10 A. J. S. Hamilton and M. Culhane 


The index 7 in equations (^) and (^) is the same 7 as in the representation (|^) or of the correlation function in 
eigenfunctions . The factor of a half which multiplies the index 7 in equations (|41^) and arises because the 

correlation function is the square of the density. The reality conditions S’‘*{r) = hence = 5‘’{—k), along with the 

usual properties = ( — )’^Yi^-rn and Yim{—k) = { — YYimlk) of the spherical harmonics, imply 


rs* _ / \m ts 5 ®* _ / \l + m?s 

^ujlm V / ^ — uj ,l ,—m ? ^ujlm \ ) ^ — u ,—m • 


(47) 


Equating the expansion of S‘’{r) in equation to the Fourier transform of the expansion of Y^k) in equation shows 
that the eigenfunctions and are the same up to factors C{'y/2+iuj,l)-. 


— C{'yl2-\-iLO, 1) S^lm 


C{r,,l) 




r[^] 
r[^] ■ 


(48) 


The r functions come from integrals of power laws with spherical Bessel functions. The two T functions in equation (^^ can 
be reduced to the product of a single T function with a sine function and a rational function, but the expression (^ ) manifests 
the symmetry between the two representations S'lid ^Zim¬ 
in wZm-space, with 


rj = 7/2 + iu> 

the spherical distortion equation (^) becomes 
(2-7)(1-?7 + q) 


1 + /? 


(2-?7)(3-??)-Z(/ + 1) 


(49) 


(50) 


with an identical relation between 5^;^ and Suiim- As expected, the spherical distortion operator reduces to eigenvalues in 
wZm-space. 


4.3 The correlation function of logarithmic spherical waves 

As will be seen shortly, the normalization of the ‘Fourier’ representation S^i^, equation ( ph), l eads to simpler expressions for 
the correlation function in a;Zm-space than that of the ‘real’ representation SZimj equation (^). We therefore particularize to 
the former. 

Statistical isotropy about the observer implies that the correlation function of logarithmic spherical waves is diagonal 
in the angular indexes Zm: 


mi ) ( ) ^^1^2 ^mi,—m 2 


(51) 


where the Jq^and Smi,-m 2 denote Kronecker deltas, not to be confused with the overdensity. The extraneous minus signs 
in equation (^^ arise from pair exchange symmetry and reality, as in equation (^), and would disappear if the correlation 
function in equation © were defined by {S°Y*) rather than {S°Y}. Equation (il) is valid for both redshifted and unredshifted 
correlation functions. 

The unredshifted correlation function from equations (0), (0), (0, and (0, 


^uJlUJ2l ( ) (,^uJil7n^uj2,l, — Tn') 




fo(fel +fe2)P(fcl)fcr^^"'“'ZC2""'/^“‘“"(-)Vim(fel)Y)m(fc2)d®fcld®fc2 . 


(52) 


Equation (0) reduces to the simple result that the unredshifted correlation function ^,., 11 . 12 * i® a function only of the sum 
Wl + CU 2 : 


e 


UJ 1 UJ 2 I 




ijJl2 = Cell + Ctl2 


(53) 


where the reduced correlation function in equation ( p3[ ) is the representation of the power spectrum P{k) in eigenfunctions 
^-3+7+i“12 


P(fe) = 




^aji2 (^Tr) 


P{k)k- 


,3-7-11^12 


dk/k . 


(54) 


Equation (0) is entirely analogous to the representation (0) of the correlation function in eigenfunctions ■ Equating 

the expansion of ^(^ 12 ) in equation (0 ) to the Fourier transform of the expansion of P{k) in equation (0) shows that the 
‘real’ ^^112 and ‘Fourier’ ^„i 2 unredshifted correlation functions differ only by a factor 

^LU12 = (271 -)"^^^^(7 + ia;i 2 ,0) iuji2 (55) 

where C is the function defined in equation (|48|). 

It follows from equations (0), (0) and (0) that the unredshifted correlation function of spherical waves in 

logarithmic real space is related to the unredshifted correlation function ^ 0^12 defined in equation (| 2 l|) by 
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_ / \7n / X X \ / ^^/'o__^3/2 5 0 1^2, 0 

“I ) \^LJilm.UuJ2,l,-m) \ ) 14, 

C (7 + 10412, U) 

which is more complicated than the corresponding simple ‘Fourier’ relation 
The equality 0112 = o4i + 012 is a consequence of the equality of operators 


(56) 


d 


9 In ri 2 


A 


d d 

31nri (llnr2 


(57) 


Equation (^^ can be interpreted as meaning that an expansion (by a factor) of the separation ri 2 at fixed triangle shape is 
equivalent to a combined expansion of the legs n and r 2 at fixed angle between them. An equivalent statement is valid in 
Fourier space. 

In o4fm-space, with 


rji = 7/2 + io4i , 72 = 7/2 + io42 (58) 

the spherical distortion equation (^^ becomes 

_ /_ \l + m /^s ?s \ 

^uJiuJ2l V / \^uJilm^uJ2,l, — Tn/ 


1+/3 


(2 - 77i)(l - rji+ a) 
(2-7i)(3-r?i)-Z(Z + l) 


1 + /3 


(2 - ? 72)(1 -rj 2 + a) 


{2 - V2){3 - V 2 ) - l{l + 1) 


■ 


In the plane-parallel limit, wi « —UJ 2 —> 00 and I 00 , and the spherical distortion equation (p9|) reduces to 


C^u.2l = (1 + 


Ujf 


Ji 


Lof + P ^2 + 

which is Kaiser’s equation (W), except that the power spectrum is expressed in the representation equation (b4[). 


(59) 


(60) 


5 MEASURING OMEGA 
5.1 Considerations of optimality 

The spherical distortion equation ( p^ in the 0112 A-representation is a linear equation in the unknown quantities 

^“(0412, A) = Ai(a4i2, A)/3 *^i.,^2 (61) 

where indices i run over 0, 1, 2, and we adopt the summation convention for these indices. The spherical distortion equation 
( pol) in the wZm-representation can be cast in a similar form: 

I"( 0412 , A) = Ai(o 4 i 2 , A)/ 1 ‘|„j 2 ( 62 ) 

where A is shorthand for the ‘shape’ variables [(o 4 i — a42)/2, Z] at fixed 0412 = tJi -I- 042 . The quantity on the left hand sides 
of equations M and ( 1^ is the redshift space correlation function, an observable quantity. On the right hand sides of 
equations (§) and ( 1^ is the ensemble average distortion predicted when the underlying pattern of clustering is statistically 
homogeneous and isotropic. The problem is to find the values of I3'^uii2 that give the best fit between observation and 
prediction. The ratios of any two of the three fitted quantities will then give values for (3 that are 

independent of convenient to think of /3 and 0^ as ‘independent’ parameters to be fitted to the data, which can be 

combined into a single best fit at the end of the day. 

For a sufficiently large survey, the maximum likelihood solution of equation (§) at any particular 0412 is the least-squares 
solution of (Kendall & Stuart 1967, section 19.17) 

X-JJ ir (W12, A) - Mu;i2, a') - A0LO12, A')P0^,,]WiLOi2, A, A')dAdA' (63) 


where the weighting function IF( 0412 , A, A') is the inverse of the covariance matrix 

lF(a4i2, A, A') = (Af (0412, A)Ar ( 0412 , A'))"' . (64) 


The volume element dA of triangle configurations is 


dA = 


ridrir 2 dr 2 


(65) 


in terms of which the pair-volume element is d®rid®r 2 = Srr^dA rf 2 dri 2 /ri 2 , the factor coming from integration over 
orientations. Similar equations are valid in the o 4 Zm-representation (just put tildes on A, A and W in equations and 
Ip^), the volume element dA of configurations being 

dA = d(a4i— 042)72 ( 66 ) 
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with integration implying summation over angular indices Im. For complex as here, the integrand of equation (^) should 
be interpreted as representing its real and imaginary parts separately, or else some combination thereof, according to the form 
of the covariance matrix in equation (164) . 

The sum of squares equation (p3|), is minimized when the derivatives with respect to the parameters are zero, 

= 0. This is a set of linear equations, whose solution is 

= M-/{uJr 2 ) J A-(wi 2, A)r (a;i 2 , A')VF(a;i 2 , A, A')dAdA' (67) 

where M^“^(tJi 2 ) is the inverse of the symmetric 3x3 matrix 

Mij{uJi2)= Aj{ui2,A)Ai{L0i2,A')W{uJi2,A,A')dAdA' . (68) 


If the covariance matrix is diagonal, then it is trivial to invert and hence determine the optimal weighting function W, 
equation (^). As made clear by Feldman, Kaiser & Peacock (1994) and exploited by Heavens & Taylor (1995), this is one of 
the great advantages of working in Fourier space for Gaussian fluctuations, that the independence of (unredshifted) Fourier 
modes causes the covariance matrix of the (unredshifted) power to be diagonal. Feldman et al.’s minimum variance pair 
weighting [1 + 'l'P(A:)]“^ is dehnitive at scales which are Gaussian but smaller than the scale of the survey. 

As argued in Section the advantages of mode independence are preserved at least partially in the a;Zm-representation. 
We pursue this idea momentarily, before abandoning it in Section 5.2 in favour of the tJi 2 A-representation. The covariance 
matrix of the correlation function in a;Zm-space is 


(AC 


(69) 


where i)® is the four-point function, which vanishes for Gaussian fluctuations, and is neglected hereafter. For the unredshifted 
correlation function, equation ([s^, the covariance matrix (^^ reduces to 




(70) 


For a power-law power spectrum, the unredshifted correlation function in oiZm-space is a Dirac delta-function, C 12 = fc(<^i 2 ). 
In this case, the only non-zero elements of the covariance matrix (^) are the diagonal elements {A^uiiuj 2 imA^^_^^^i^), and 
these elements are all equal. In effect, the covariance matrix is the identity matrix. The covariance matrix of the redshifted 
correlation function is related to its unredshifted counterpart by the spherical distortion equation (^). In redshift space, the 
covariance matrix (^) remains diagonal for a power-law power spectrum, but it depends on /3. 

The above argument, that the covariance matrix ( p9|) is almost diagonal in the ojZm-representation, contains a serious 
flaw. Namely, it neglects to take into account the shot noise caused by the discrete sampling of galaxies. Shot noise destroys 
scale invariance, and introduces off-diagonal correlations between modes. We are not sure whether this difficulty is fatal. 
Whatever the case, it is enough to persuade us to switch to the 0112 A-representation, which proves more tractable. 


5.2 0112 A-representation 


Rather than attempt a weighting that is absolutely optimal, let us instead adopt a weighting function W that is diagonal in 
real space. Great precision in the choice of W is not required, since any weighting in the vicinity of minimum variance should 
give a variance not much different from the minimum. Thus we seek to minimize at each separation ri 2 the sum of squares 


^(’' 12 , A)-/ Ai(oii 2 , A)/3‘^„42r( 


^du!l2 


W(ri2, A)dA 


(71) 


A near-minimum variance form of the weighting function W (ri 2 . A) is (Hamilton 1993b, section 5; notice the following formula 
is a volume weighting not a number weighting, which accounts for the 4>(ri)4>(r2) in the numerator) 


IF(ri 2 ,A) 


_■I>(ri)<I’(r2)_ 

[1 -I- 'l>(ri) J(r 12 )] [1 -I- 'l'(r 2 )J(ri 2 )] 


(72) 


where J{ri 2 ) ~ ^dK. 

According to equations ( p^ and (0) , the coefficients Ai(aii 2 , A) separate into sums of products of the five shape-functions 
Bj{A), equations (p8|)-(p^, with functions 0 ^ 1 ( 0112 ) which depend on 0)12 but not on triangle shape A: 


Ai(a;i2, A) = Bj{A)aji{uJi2) ■ 


(73) 


It is convenient to imagine the five quantities as ‘independent’ parameters to be htted to 

the data. Minimizing the sum of squares ([nl) with respect to these parameters yields hve equations (cf. equation |^): 


Mi/(ri2) 


Bi(A)^'’(ri 2 , A)W(ri 2 , A)dA = / aii(a)i2)/3*C‘^i2f’i 


^dLUl2 


(74) 


where M^^^(ri 2 ) is the inverse of the symmetric 5x5 matrix 
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Mi,(ri2) = J Bi{A)Bj{A)W{ri2,A)dA . 

Define H|(ri 2 ) to be the vector on the left hand side of equation ([74[): 
HKri2) = M^i(ri2) / i3,'(A)r (ri2, A)W(ri2, A)dA . 


(75) 


(76) 


As will become apparent below, the H|(ri 2 ) can be interpreted as the generalization to spherical distortions of the harmonics 
of the redshift correlation function in the plane-parallel case. On the right hand side of equation ([t^, the quantities aji{uji 2 ) 
can be replaced by an operator in real space and taken outside the integral. Equation (M) then becomes, explicitly. 


Ho(ri2) = (^1 -h |/3 + i/3^) ^(>" 12 ) 
S|(ri 2 ) = + y/3^) [C(fi 2 ) - ^(’' 12 )] 


HKn 2 ) = 


C(ri 2 ) -h ^C(f’i 2 ) - ^?(ri 2 ) 


=^^■ 12 ) = i/3aC(’'i2) +13^ |iaC(f’i2) 


-h a 




?(’' 12 )] } 


H3(ri2) =/3^||^(ri2)-^C(»'i2)-a ic(7-i2) + :j^C(ri2) - [^I(n2) - ic(n2) - ^^(ri2) | 

where, following Hamilton’s (1992) notation, 


(77) 

(78) 

(79) 

(80) 
(81) 


|(ri 2 ) = 


(82) 


2 /*^12 _ r'f'12 

f(r)rdr , f(ri 2 ) = / ^{r)r‘^dr , f(ri 2 ) = 5rC/ / 5(r)r‘‘dr . 

Jo Jo Jo 

In the plane-parallel limit, Hq, and H| go over to the monopole, quadrupole and hexadecapole harmonics of the 
correlation function, while Hf and vanish. The three equations ([^-([^ for Hq, H 2 and HI look exactly like their plane- 
parallel counterparts, equations (6)-(8) of Hamilton (1992). In particular, the ‘quadrupole’ to ‘monopole’ ratio 


H|(ri 2 ) 


1/3+1/?" 


Hg(ri2) - .=.o(ri2) 1-1-|/3-I-I/32 


(83) 


provides a way to measure /3 in a manner independent of the shape of the correlation function, but now for fully spherical 
distortions. 

As remarked by Cole et al. (1994, appendix B), the combinations of ^ on the right hand sides of equations (^^-(^^ can 
be regarded as arising from windowing the power spectrum with spherical Bessel functions ji{kr 12 ). Thus if Hf(fc) for I — 0, 
2, 4 are defined by 


Hf (fc) = 47ri 


2 /(ri 2 )ji(A:ri 2 )ri 2 dri 2 


H/(ri 2 ) = 


27r2i* 


'^i{k)ji{kri2)k^dk 


then the three equations (f^-d^) become 


(84) 

(85) 

( 86 ) 

(87) 

an also be written 

( 88 ) 

The presence of the inverse 5x5 matrix Mr^{r 12 ) in equation @ makes the quantities (ri 2 ) more difficult to evaluate 
than in the plane-parallel limit. However, the 5x5 elements of the matrix, equation ([75|), are integrals over known quantities 
Bi{A) and the chosen weighting function IT(ri 2 , A), so can be precomputed. 


H?(fc) 

= (i + 2/>+l(j")p(/.') 

H|(fc) 


i^fc) 

= |/3"P(/c) 

where 

P{k) is the unredshifted 

ii(fe) 

|/3 + |/3' 

Hg(fc) 

l + i/3+|/3^ ■ 


5.3 The validity of the constant a approximation 

Equations ([77|)-(^l|) are valid provided that a(r), equation (j^), is approximated as a constant (cf. Section |3.l|) . How can the 
validity of this approximation be checked? 

Mathematically, the passage from equation ([hI) via equation (|^) to equations (|7^)-(^) involves the approximation that 
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j a(ri)Bi(A)B, (A)W(ri2,A)dA = aMi, (ri2) (89) 

where a on the right hand side is required to be a constant independent of ri 2 and ij. A similar constraint on a(ri)a(r 2 ) is 
required. Note that a{r) in equations ( p^ ) and multiplies only the shape functions Bi and B^, so it suffices that equation 
( p9| ) be valid when either of i or j is 1 or 3. The validity of equations should be checked in applying equations ©-(il)- 


5.4 Deconvolving the catalogue window 

In measuring the redshift space correlation function ^'’(ri 2 , A), it is essential to disentangle the true anisotropy of the redshift 
correlation function from the artificial anisotropy induced by the angular and radial selection functions of the survey. This 
deconvolution is easy in real space, where the observed galaxy density is the product of the true density and the selection 
function of the survey, but is complicated in any other space. 

The ‘deconvolution’ procedure is described in some detail by Hamilton (1993b, section 6) for the plane-parallel case where 
the aim is to measure the redshift correlation function ^'’(ri 2 ,p) as a function of separation ri 2 and cosine angle ^ to the line 
of sight. The procedure generalizes easily to the present case where the redshift correlation function ^'*(ri 2 , A) is a function of 
separation ri 2 and triangle shape A. Briefly, each observed galaxy pair at separation ri 2 and triangle configuration A makes 
a contribution to the correlation function ^'*(ri 2 . A). To correct for the selection function of the survey, one simply divides by 
the probability of finding a pair at ri 2 and A, given the boundaries and selection function of the survey. 


6 CONCLUSIONS 


The aim of this paper has been to present the theoretical foundation of a procedure for measuring the linear growth rate 
/3, hence the cosmological density H in unbiased standard cosmology, from spherical redshift space distortions in a manner 
independent of the shape of the power spectrum. 

Our most important point, presented in Section 3, is that there exists an operator d/91nri2|A which both commutes 
with the spherical distortion operator, and defines a characteristic scale of separation ri 2 . Ratios of the amplitudes of the 
eigenfunctions of this operator provide measures of f3 independent of the shape of the power spectrum. 

In Section 4, we presented a complete set of commuting operators for the spherical distortion operator. The eigenfunctions 
of this complete set are spherical waves about the observer, with radial part lying in logarithmic real or Fourier space. 

In Section 5, we discussed the practical measurement of f3 using the ideas presented. In particular, we showed that there 
is a set of five functions Hf(ri 2 ), equations ([7^-(^l|), which can be regarded as the generalization to fully spherical distortions 
of the monopole, quadrupole, and hexadecapole harmonics of the correlation function in the plane-parallel case. 

A drawback of the method is that the spherical distortion operator commutes with d/(91nri2|A, and likewise with the 
complete set of operators in Section 4, only to the extent that the logarithmic slope of the radial selection function can be 
approximated by a co nstant. We argued in Section I.l that this may be a reasonable approximation in practice, and we 
showed in Section p.3| how to check the validity of the approximation. 
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